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ENGINE ANGULAR MOMENTUM ON LONGITUDINAL AND DIRECTIONAL 
STABILITY IN STEADY ROLLING MANEUVERS 
By Ordway B. Gates, Jr., and C. H. Woodling 


SUMMARY 


The effect of engine momentiim on the longitudinal and directional 
stability of aircraft in steady rolling maneuvers has been investigated. 
The results presented indicate that the gyroscopic moments produced on 
the aircraft by rotating engine in rolling maneuvers can have an 
appreciable effect on the range of rolling velocities for which longi- 
tudinal or directional instability might occur. 


INTRODUCTION 



t 


The analysis presented in reference 1 of the effect of steady rolling 
on the longitudinal and directional stability of aircraft was made for the 
assumption of zero engine momentum^ hence, the results presented were 
Independent of the direction of rolling. Some of the present-day aircraft 
have exhibited different characteristics in left and right rolls which can 
be attributed to the asymmetric moments produced on the aircraft by the 
rotating engine. The purpose of this report is to present the aircraft 
equations which include these asymmetric engine gyroscopic moments and to 
demonstrate the effects of these terms on the divergence bo\mdaries pre- 
sented in reference 1 for the steady rolling case. The divergence bound- 
aries presented in this paper are for aircraft having static stability. 


SYMBOLS 


W weight, lb 

T total aerodynamic moment, Ib-ft 


I 
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L rolling moment, lb -ft 

M pitching moment, lb -ft 

N yawing moment, lb -ft 

H angular momentum, slug-ft^/sec 

moment of inertia about body X-axis, slug-ft^ 

p 

ly moment of inertia about body Y-axis, slug-ft ■ 

p 

I 2 moment of inertia about body Z-axls, slug-ft 

lyz product of inertia in X,Z-plane (positive when principal 

X-axis is below body X-axis at nose), slug-ft2 

Iv moment of inertia of engine about body X-axis, slug-ft^ 

e 

p rotational velocity about body X-axis, radians/sec 

q rotational velocity about body Y-axis, radlans/sec 

r rotational velocity about body Z-axls, radlans/sec 

03g engine rotational velocity, radians/sec 

V aircraft velocity, ft/sec 

u component of V along X-body axis, ft/sec 

V component of V along Y-body axis, ft/sec 

w con 5 )onent of V aAong Z-body axis, ft/sec 

a angle of attack of X-body axis, w/u, radians 

^ angle of sideslip, v/u, radians 

F total aerodynamic force, lb 

Fx. Fy, F^ Components of the aerodynamic force along X-, Y-, and 
Z-axes, lb 

Zj, mj, nj. Direction cosines relating aircraft body axes to earth 
Z-axis 
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P 

^Q’ 

aij., a.-^, ^2) aQ 



damping ratios in pitch and yaw, respectively 
coefficients of characteristic equation 


Dt 


total time derivative 


A dot over a symbol .indicates differentiation with respect to time. 


EQUATIONS OF MOTION 


The aircraft equations of motion written in vector form are: 


^(H) = T 
Dt 


mJ5_(v) = F + W 
Dt 


( 1 ) 


k 
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where — refers to total differentiation with respect to time, H is 
Dt _ 

the angular momentum vector, T represents the aerodynamically applied 
moments, V is the aircraft velocity vector, F represents the aerody- 
namically applied forces, and W is the weight of the aircraft. A 
right-handed system of axes is chosen which originates at the center of 
gravity of the aircraft and which is fixed in the aircraft. The X-axis 
is assumed to he coincident with the X-axis of the engine and the 
X,Z-plane is considered to he the plane of symmetry of the aircraft. 
Also, the mass dlstrihution of the engine is assumed to he symmetrical 
about the X-axis, and the engine is assumed to he rotating with constant 
speed. The rotational velocities about the X-, Y-, and Z-axes are p, 
q, and r, respectively, and the components of V in this system of 
axes are u, v, and w. 

For the previous conditions, the vector momentum is given by 


H 


= K- 


ix? 


■ ^XZ^ 






+ + ki 




where i, j, and k are unit vectors in X-, Y-, and Z-directlonsj I^, 
ly, and are the moments of inertia; is the product of inertia 

in the XZ-planej is the moment of inertia of the engine about the 

X-axis; and o)g is the rotational velocity of the engine about this 
axis, taken positive in the same sense as the rolling velocity p. 

Equations (l) become, after differentiation and resolution into 
X-, Y-, and Z- component s : 

Rolling : 


" ^xz^ " (% “ - ixzPG ~ 


(2a) 


Pitching: 


lyG 


iz - 


lylpr + 1 


XZ 


(p= - 






M 


(2b) 
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Yawing: 







N 


(2c) 


X- force : 


m( u + qw - vr ) = ^ + 


WZ. 


(2d) 


Y-force; 


.(v - pw + ur) = ^ Fy + Wmj 


(2e) 


Z-force: 


m(w + pv - uq) = \ F^ + Wn^ 


(2f) 


The terns Zj, m^, and n^ are the direction cosines between the earth 

Z-axis and the axes being used, which are fixed in the body. The equa- 
tions which relate m^, and n^ to the airplane rotational velocities 

p, q, and r are 

Zj = m^r - n^q 
m, = n^p - Z^r 

5 5 5 

^3 = ^3^ - 1 Q 3 P 

Certain assumptions beyond that of a constant rolling velocity are 
necessary in order to linearize these equations. The term p2 - r2 in 
the pitching equation is taken as approximately equal to p^, the term 
I;j^qr in the yawing equation is considered negligible, and it is 
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assumed that no change occurs in the X-component of the forward velocity. 
Further, it is assumed that 




a 


and, for the assumption of small out-of-trim aerodynamic forces, the 
Y-force and Z-force equations are approximately 

P = pa - r 


d = q - pP 

Also, the aerodynamic moments are taJien as 

M = ^ Zee = M q + Mp, Ax 

dq oa liu- 

N = + ^p = Nj-r +■ KeP 

Br 3p 


Equations (2), for these assumptions and substitutions, become in determi- 
nant form 


M. 


D 


-(^X - %)Po ■ ^Xe"e 


Ax 


P 

0 

Np 

0 --E- 


Iy 


(ix - Iz)pq + 


Iv 


- ^XZ 2 

T 

-*-Y 


= 0 


0 

-1 


-P. 


D 


1 

0 


Po<^o 
= 0 


where p^ refers to a constant value of the rolling velocity p and 

ttQ, to the initial value of the angle of attack. Expansion of this 

determinant with the right-hand side set equal to zero yields a charac- 
teristic equation of the form 
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4- + Sq ~ 


(3) 


where 


aij. 


1 


B.X = 


^ ^ 


a = t)2+!!P_^ + 

2 ° Iz % %^Z 


(^Z " ^x)Po ■ 


(% " %)Po + ^Xg'^e 


= -Hr p 2 ^ 2 ^ 

IZ ° ly ° lylz lylz 


_ p N 

lylz I; 


(iz - Ix)Pq - 


lylz ° ly 


(% - %)Po + %e“e 


(^Z “ ^x)^o ^Xg*^e 


(% “ %)Po + ^Xg'^e 


For the substitutions 




2 _ 


Vo^ 


U), 


2 _ ~^cx, 


IvP “= 
y-^o 


-M. 


2| otUo - 


■00 Ip 
•^y-^o 
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= 


-Nr 

^Zf>o 


T 




P 


O 


These coefficients hecome 
= 1 

aj = Po(2^y^^ + 2|g0ig) 


ag - 


P P - % - %e'^)Px - ly + ly T 

1 + oi^2 ^ ^ _e_M e^ + 


)(^ 






'1 = Po^(2^“> + 2|gaig + 2|gCSga>4,2 + 2|^C0^g^ 


ao - Pc 




I 7 - Iv - 


2 2 

+ a>, 0)„ + 

i|/ 0 



which are essentially equivalent to the coefficients presented on page 9 
of reference 1, with the exception of the engine momentum terms. It will 
he noted that the Pq factors which multiply the coefficients shown 
here do not appear in reference 1. This is attrihutahle to the fact that 
the differential operator used in this paper is a time operator, whereas 
in reference 1 the operator has been made a function of rolling velocity. 
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ANALYSIS AND DISCUSSION 


The combinations of co, and which will res\xlt in the constant 

\|r 0 

term of equation 3 a_ being negative, for given values of i., 

U u Y 

and T, and hence give aperiodic instability are of primary interest. 

For an airplane with given values of pitch and yaw natural frequencies, 

these values of go, and go„ define the rolling velocities for which 

y U 

this instability will exist. 


A sample of these aperiodic divergence boimdaries is shown in 
figGire 1 for t = 0(aOg =0), for = 0 and = O.OO 3 I. The mass 

and aerodynamic characteristics of the airplane for which these bound- 
aries were constructed are presented in table I. It should be noted 
that boundaries constructed for a constant value of , do not corre- 


spond to a constant value of 


M N^ 

— — , but, rather, every point on the 


curve represents, dimensionally, a different value of this parameter. 

In this plane, the pitch and yaw frequencies of a given airplane, for 
all values of lie along a straight line similar to the one shown 

in the figure. The point shown for p^^ = 1 radian/sec defines the fre- 
quencies of the airplane chosen for this illustration. The slope of 
this line is determined from the ratio of the square of the pitch and 
yaw natural frequencies. For the case shown, the frequency locus of 
this airplane passes through the divergence boundary constructed for 
= 0.0031 for Pq = ±1.8 radlans/sec and remains on the unstable 

side of the boundary up to p^ = ±2.3 radlans/sec. Generally, the char- 


acteristic roots of the system in the unstable region of this plane are 
a pair of stable complex roots, one stable real root, and one unstable 
real root. 


It would be possible to take into account the engine moment-um by 
plotting boundaries for various values of but, as was mentioned 


for the I , f 
> 0 


case, these boundaries would not correspond to a constant 


value of engine momentum. Also, both positive and negative values of 
I^ T would have to be considered in order to cover both the left and 

right rolling conditions. 


10 


NACA RM L55G05 


The former difficulty can he avoided for both these cases by 
plotting the boundaries as a fimction of the dimensional frequency param- 
eters rather than in terms of a>^,2 and a)g2. This necessitates the Con- 
struction of a boundary for each rolling velocity^ but this construction 
is relatively simple. A sample of these boundaries is shown in figure 2 
MnNr* , » 

for — — = 0.044 and ly ^e = 0* For this case the boundaries for left 

and right rolling are identical. The effect of the engine momentum on 
these boundaries can be seen in figure 5 for tiio case of 

Ix^cOe = 17^55^ slug-ft^/sec. Boundaries are presented for both right and 

left rolls^ and the critical roll velocities are shown on the figure for 
both rolling conditions. For positive (right) rolling the unstable range 
of p is between p^ = 2.1 radlans/sec and p =2.5 radians/sec and 

for negative (left) rolling this range is between p^ = -1.7 radians/sec 

and p^ = -2.2 radlans/sec. In order to specify the absolute range of the 

roll rate vhich might be critical for a given airplane, it is necessary to 
know the critical values of p for both the left and right rolls. For 
the value of Ijr considered here, this unstable range for the example 

airplane would be defined as being 


1-7 < |Po| < 2-5 


Curves are plotted in figure 4 which show the effect of engine momentum 
on this absolute range of critical p for the particular airplane being 
considered. For oUg = 0, this band is between jpj of 1.8 and 

2.3 radians per second, as mentioned previously, and broadens as the 
engine momentum increases. Thus, the range of rolling velocities for 
which a given airplane might ejtperlence instability, based on this steady 
rolling assung)tlon, can Increase appreciably with the magnitude of the 
momentum of the rotating enginej hence, the effect of engine momentum 
should be considered in the analysis. 


A point of interest with respect to the construction of the diver- 
gence boundaries in the dimensional frequency plane is that it is possible 
to obtain mathematically the envelope of these boundaries by plotting the 
locus of the points of maximum cunvature of the respective curves. The 
mathematical expression for the curvature of a given function can be found 
in any calculus book (for example, ref. 2) and in order to obtain the 
desired envelope it is necessary only to maximize this expression in the 
proper manner. From this envelope and the zero and asymptotic values of 


W, 




and 




which are rather easy to calculate, one can approximately 



MCA RM L55G05 


11 


draw the divergence boundaries for any given value of roll rate for 
specified inertia characteristics. In any event, the envelopes of the 
branches of the curves will define, for the steady rolling case, the 
combinations of pitch and yaw frequency for which there will be no roll- 
induced instability. The equations from which the envelopes can be 
calculated, and expressions for the previously mentioned asymptotic values 


Nc 


of 


and 




are : 


Np 


(a) Equations for determination of zero and asymptotic values of 


and — : 




Po%3 


(^Z ' ^x)Pc 



M(^Nr 


[y 2 ^Xg^ePc 

-Po + 


+ 


IXg^ePo 


"X 


iz 



1 

11 

(ix - Iy)po^ + IXg^ePo 




ly 


/ . _ 2 

~ lylz Ax - WV 2 . ^Xe">ePo 

' _ (^Z ~ ^x)pq^ V % / H 

ly ly ly 
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(l) Equations for determination of envelope of divergence boundaries 






+ 2 ^ 




p = a ± b 


Ma _ 


IXg^ePo - (^Z - Ix)Po^ 




/M W 

l^iyiz ° 


The equations presented for determination of the envelopes of the diver- 
gence boixndaries require some further explanation. The combinations of 
Np -Ma 

— and which define the envelopes are 


and 



a + b 



- b 



a - b 



c + b 


The branch of the family of divergence boundaries to which each combina- 
tion applies depends on the sign of the rolling velocity. A sample of 
the envelopes calculated for the boundaries of figures 2 and 5 is shown 
in figure 5- 
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CONCLUDING REMARKS 


From the analysis presented in this paper it appears that the 
momentum of the rotating engine of an aircraft can increase appreciably 
the ranges of rolling velocity for which the aircraft might experience 
a roll-induced aperiodic divergence in steady rolling maneuvers. For 
the aircraft discussed herein^ the range of critical rolling veloci- 
ties p was calculated to be 


1.8 < p < 2.5 


when engine momentum was assumed zero, and 


1-7 < 



< 2.5 


when engine momentum was considered. For higher values of engine momentum 
than that assumed for this illustrative example, this band would, of 
course, be further expanded. 

It was also shown how the construction of the divergence boundaries 
was affected by inclusion of engine momentum, and an alternate method of 
construction to that presented in NACA TN 1627 was discussed. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., June 50^ 1955 • 
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TABLE I.- MASS AED AERODYNAMIC CHARACTERISTICS 
OF EXAMPLE AIRPLANE 


Ijj, slug-ft^ . . . 

ly^ SlUg-ft^ . . . 

I^^ slug-ft^ . . . 

P 

Iy2, Slug-ft . . 

ly slug-ft^/sec 

% ^ 

sec 

ly sec 

Mg 1 

% sec^ 

Np 1 

sec^ 

V, ft/sec . . . . , 


10,976 

57,100 

64,975 

942 

17,554 


- 0.105 


- 0.421 


-5.30 


2.38 

691 
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1.- Boundaries in the ci:q^ plane which define regions of 

aperiodic divergence for example aircraft. 
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-Mg 1 

I Y sec 2 


Pq, rad ions /sec 


1.0 1.5 2.1 2.5 



I 

’ sec2 


(a ) Right rolls . 


Figure 5*- Boundaries in the Hx/^Y plane for exan^jle aircraft 

which define regions of aperiodic divergence as a function of rolling 

.2 


velocity. = 17 ^ 55 ^ 


slng-ft 


sec 


18 


Pq, radians /sec 


unstable 


(b ) Left 


Figure 3-- Coi 
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Figure 5 .- Envelopes of divergence boundaries presented in figures 2 and 5- 
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